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Abstract 

In this paper we discuss the systematics of quarkonium production at the LHC. In particular, 
we focus on the necessity to sum logs of the form log(Q/pj_) and \og(px/mn). We show that the 
former contributions are power suppressed, while the latter, whose contribution in fragmentation 
is well known, also arise in the short distance (i.e., non-fragmentation) production mechanisms. 
Though these contributions are suppressed by powers of mg/pj_, they can be enhanced by inverse 
powers of v, the relative velocity between heavy quarks in the quarkonium. In the limit p± 3> mg 
short distance production can be thought of as the fragmentation of a pair of partons (i.e., the 
heavy quark and anti-quark) into the final state quarkonium. We derive an all order factorization 
theorem for this process in terms of double parton fragmentation functions (DPFF) and calculate 
the one- loop anomalous dimension matrix for the DPFF. 
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Quarkonium production is a semi-inclusive hadronic observable that requires minimal 
non-perturbative input. Predictions for this observable depend only on the usual parton 
distributions and a set of local non-perturbative quarkonium production matrix elements 
that can be extracted from the data. These predictions are based on the Non-Relativistic 
QCD (NRQCD) [l| factorization theorems and are formulated as a double expansion in a s 
and v, where v is the typical relative velocity of the heavy quarks in the bound state. 

In this paper we will concentrate on the vector states J/ip and T, which have the quan- 
tum numbers 2S+1 Lj = 3 S\. At leading order in the v expansion there is only one relevant 
NRQCD matrix element, which represents the probability of a heavy quark-antiquark pair 
in a color- singlet 3 Si state to form a quarkonium bound state. However, sub-leading con- 
tributions in the velocity expansion can receive kinematical enhancements [2] that scale as 
powers of p±/rriQ. Thus calculation of high p± quarkonium production should be formulated 
as a systematic expansion in three parameters: a s ,v, and m,Q/p±. The relative importance 
of the various mechanisms depends on all three parameters. The parameters a s {2mo) and v 
are fixed for a particular quarkonium state, but p±/rriQ clearly varies within the experiment 
depending on what p± is measured. Currently, the LHC experiments CMS and ATLAS have 
measured J / ip production in pp collisions at _Js = 7 TeV with p± as high as 70 GeV (3), Q 



The dominant contribution at asymptotically large p± will come from gluon fragmenta- 
tion, shown in Fig. [T]d. In this paper we will follow the terminology of Ref. {tJ and refer to 
direct production via non-fragmentation processes as short distance (SD) production. As 
will be explained below the fragmentation contribution naively scales as (pj_/ m o)' y4 relative 
to SD production contributions shown in Fig. [T^i. The v 4 suppression is due to the fact that 
in gluon fragmentation the quark pair is produced in a color-octet state and color quenching 
requires subsequent emission of soft gluons which vanishes in the static limit. But, as will be 
shown below, this SD contribution is actually further suppressed by a factor of (mg/pjj at 
leading order in a s (p±). Thus the fragmentation contribution scales as {p\/^q)v 4 relative 
to LO SD production. However, the NLO color-singlet SD contribution scales as and so 
fragmentation scales as (p 2 ± / 'mg)(v 4 / 'a s ) relative to it. While this contribution is suppressed 
by an additional power of a s (p±) relative to the LO SD color-singlet piece, it does not fall 
off as steeply with p±. Thus at lower p± ~ 2mg we expect the color-singlet SD mechanism 
to dominate, while at very high p± ^> 2mQ we expect fragmentation to dominate. The 



and T production with p± up to 24 GeV 
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important question, however, is how each contribution affects the differential cross section 
in the intermediate p± region. 

A crucial distinction between high and low p± predictions is the size of the large logs that 
arise in perturbation theory. When the hierarchy 

Q>p ± >m Q (1) 

exists, both log(Q/p±) and log(p^/mg) may appear in perturbation theory, where Q is the 
underlying hard scattering scale. The logs of p±/rriQ in single-parton fragmentation can be 
resummed using standard renormalization group (RG) techniques. This has already been 
accomplished in the literature [t 9 ] . However, the same type of logs arising in SD production 
have yet to be summed. These logs can be resummed by thinking of SD production as arising 



from double parton fragmentation (DPF) and then utilizing standard RG techniques 10 ]. 
Of course, it is important to remember that the growth of these logs is accompanied by a 
power suppression and thus we expect these logs to be numerically important only in the 
intermediate, as opposed to asymptotic, regime. Phenomenologically, most of data on J ftp 
and all data on T production at the Tevatron is at moderate p± where both single parton 
fragmentation and SD production mechanisms are important. This will continue to be true 
at the LHC, except for the very highest p± for J/ip production, where one is plausibly in 
the fragmentation regime. 

Let us determine the size of the logs in the intermediate regime where most of the available 
data is. Fragmentation and SD production are of the same order when 1 

4 2 4 

P± v 4 ~i (LO) , o P± . 2 — ~ 1 (NLO). (2) 



(2m Q ) 4 v ' (2m Q ) 2 a 

For the J/ip (T) system we will take a s ~ v 2 ~ 0.3 (0.1). Then both expressions in Eq. (J2J) 
yield p±_ ~ 5 (30) GeV. For these values of p±, \og(pj_/mQ) is not a huge logarithm. However, 
for both Jjij) and T the ratio p±/m,Q is comparable to the ratio m^/AgcD- Resummation of 
log(m&/AQCD) is required for accurate prediction in many processes involving heavy quarks. 
It is likely such a resummation will be useful for quarkonium production as well. 



1 This is correct for comparing color-octet fragmentation to SD production. For 3 Si color-singlet fragmen- 
tation, v 4 should be replaced with a 2 but since numerically v 2 ps a s the estimate for p±/rriQ holds for 
this case as well. 
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An accurate prediction for this intermediate regime is important since the norm of the 
production cross section is fixed by the NRQCD matrix elements, which have to be extracted 



from the data. These matrix elements have 
B factories ll|, e + e~ annihilation at LEP 



jeen extracted from production of JMat the 



12| . photoproduction [l3|, B decays [l4j, Il5| . 
as well as from fixed target [16| and hadroproduction Q, Q. For a recent analysis that 



extracts NRQCD matrix elements from a g 



NLO theoretical calculations, see Refs. 



19 



obal fit to a wide range of ex per iments using 



20] . The extraction in Refs. [171 . Il8| involves 



interpolating between (non-resummed) SD production and fragmentation and it is unclear 
how to estimate the errors involved in this process given the merger regime is contaminated 
by the aforementioned large logs. Thus the accuracy of the extraction will be enhanced by 
the log resummation studied here. As such, this paper can be thought of as a continuation 
of the study started in [l?], [l8] at leading logarithmic order. 



In the asymptotic regime where 



4 

P± 



V>1, 



fragmentation dominates. A classic prediction of NRQCD in this regime is that quarkonium 
production is dominated by color-octet fragmentation and the quarkonium is produced with 



purel y tr ansverse polarization 



2lM23j|. Presently, there is no indication of this trend in the 



data 24] . Currently available data does not probe the asymptotic regime for the T system, 



but the failure of the prediction for J/ip calls into question the validity of the NRQCD power 
counting in charmonium. In the J / ip one must recall that treating this system as Coulombic 



is questionable and might require a power counting distinct from NRQCD [25|, |26[. Since 
much of the data falls in the intermediate regime we expect the resummation of \og(p±/mQ) 
to shed some light on the polarization puzzle. 

The purpose of this paper is to derive factorization theorems and evolution equations 
that will make it possible to resum log(p±/rriQ). We will also show that contributions with 
\og(Q/p±), where Q 3> p±, are power suppressed. Resummation of logfjo^/mg) requires the 
introduction of the power suppressed double parton fragmentation function (DPFF), which 



was first introduced in Ref. 



10] • We will show that the former types of logs are suppressed by 



powers of p±/ Q. We will use soft-collinear effective theory (SCET) 27], [28] for our derivation 
of factorization. 2 In Appendix A we review the minimal SCET formalism and notation used 



2 Like all SCET based factorization, our result relies upon standard methods of proof [3] when it 



comes 




b 



FIG. 1: a) The short distance production mechanism with the quark pair produced in a color- 
singlet state, b) Gluon fragmentation production contribution to hadro-production. The two soft 
gluons emitted via El transitions are suppressed by v 4 . 

in this paper. We do not explicitly perform the resummations in this paper, leaving that for 
a subsequent publication. The new results in this paper include the necessary factorization 
theorems for both types of resummations as well the one-loop anomalous dimension matrix 
for the DPFF. 

I. FACTORIZATION 

The inclusive differential cross section for the production of a quarkonium state H with 
mass Mh, four- momentum p and transverse momentum p± via the collision of two incoming 
hadrons, hi and h 2 , with momentum pi and P2, respectively, is: 

^(h l + h 2 ^H + X) = ] , 2 tEE / TT^i^iPi+P^-P-Px) 
dp{ 2( Pl + p 2 f 4 ^ ^ J (4tt) 2 

x\M(hh 2 ^ H(y,p ± )+X)\ 2 . (4) 

Here the azimuthal angle of H is integrated over, and y is the rapidity of H, which is 
restricted to be in the range —2.4 < y < 2.4 in the LHC experiments jsl, Q]. The hadronic 
matrix element above receives contributions from all scales between Q, the invariant mass 
of the partonic collision, and the hadronic scale A ~ 1 GeV. In addition to these scales there 
are two other relevant scales that may be hierarchically separated from Q and A, namely 
p± and the heavy quark mass nig. In principle, perturbation theory could be plagued by 
large logs of the ratios of these scales. Thus to be able to calculate within a well-defined 

to the cancellation of the so-called Glauber contributions. 
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approximation scheme, we need specify which regime we are in. We need to consider the 
two possible hierarchies 

(I) Q > p ± > m Q , (5) 

and 

(II) g~p±>m Q . (6) 

We will now show that the regime (I) is power suppressed. In this region, the physical 
picture is that the quarkonium must be accompanied by at least two nearly back-to-back 
jets whose net p± <C Q, and whose total invariant mass is ~ Q. Let us consider the 
appropriate theory below the scale Q. In the lab frame the incoming hadrons move along 
the z axis and their light-cone momenta (k + , k~, k±) scale as (y/s, 0, 0) and (0, y/s, 0). Since 
the quarkonium's p± is the infrared scale, the initial-state radiation collinear to the incoming 
beams has four-momenta scaling as 

Tp 1 Tp 
Pn~(Q,-^,P±) Pn~ (-q,Q,P±) • (7) 

In addition there exists soft radiation whose momentum scales as 

Ps ~ (P±,P±,P±)- (8) 

Taking as our expansion parameter A = p±/Q we see that while the collinear momenta 
scale as Q(l, A 2 , A) the soft momenta scale as Q(X, A, A). This is the scaling associated with 
SCETn, since soft modes have large enough momenta to change the p± components of the 
collinear modes, as opposed to SCETi, where the soft (in this context often called ultrasoft) 
scale as (A 2 , A 2 , A 2 ). Thus at the scale Q the hard central jets are integrated out, and we 
match onto SCETn where the infrared scale is p± and the quark mass is irrelevant (unless 
we are interested in corrections of order tjiq/Q, which we will ignore). 

Now that we have established the problem is posed in SCET n , the proof of power sup- 



pression is identical to the proof given in Ref. [30] for the case of Higgs production at 
p± rriHiggs- The suppression of central jets arises for both single and double parton frag- 
mentation, as the argument runs the same way in both cases: the two (or more) hard partons 
corresponding to the central jets cross the cut and are integrated out at the high scale as 
depicted in Fig. [2j The power suppression of such a contribution can be seen by noting that 
when there are no central jets the leading contribution scales as l/p±- This scaling arises 
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FIG. 2: A typical contribution to gluon fragmentation with two hard central jets. The lines going 
across the cut are integrated out, while the uncut lines represent the fragmenting gluon. 



either due to a collinear emission or virtual emission, which scales as S 2 (p±). When all of the 
momentum crossing the cut is hard, such a scaling factor is necessarily absent and thus the 
resulting operator which is generated is necessarily power suppressed. The situation with 
both central jets and radiation down the beam pipe (though still satisfying p± ^> ttiq) is 
reproduced by the effective theory via the one loop matrix element of the power suppressed 
operator generated by the central jets. 

Let us now consider the regime (II). When p± ~ Q the scale that controls the IR physics 
is rriQ and the scaling parameter is A = rriq / p±_ . Furthermore, we wish to keep the full 



dependence on the quark mass so we match onto SCET m [31|, |32| at the scale Q. In 
this particular case, as we shall see, there is no contribution from soft radiation at all, so 
the distinction between an SCETj and SCETn type theory becomes moot. The light-cone 
components of the quarkonium momentum scale as 



so 



P~P± 



P 




n± ) ~ p±(e y ,e y ,n±) 



with 



n 



tanh y) 



(9) 



(10) 



2 ' cosh y ' 

where n'^ 1 is given in standard four- component notation, n' ■ n' — 2, n' • p = 2m^cosh?/, 
m ± = P 2 ± + Mft, and n± is a unit vector in the _L direction. We distinguish the direction 
perpendicular to the beam direction, denoted by _L, from the direction perpendicular to the 
quarkonium direction, denoted by _L'. Corrections to these leading terms are suppressed 
by powers of rriQ/p±. The invariant mass of the final state remnants X in the production 
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process is p\ = [(pi + p 2 ) - p} 2 = AE cm (E cm - m± coshy) + M 2 H m 4E cm (E cm - p± coshy). 
We restrict ourselves to the regime where E cm > p± ~ so that the final state 

remnant can be integrated out at the hard scattering scale Q ~ p±. 3 

Since the remnants of the production process have an invariant mass p\ ~ Q 2 , they 
must be integrated out when matching onto SCET m . This is done in the manner of an 
operator product expansion where the differential cross section in Eq. fl4]) is expanded in 
terms of all the SCET operators allowed by the symmetries of the theory. We must consider 
contributions from SCET operators with light quarks and gluons that can produce heavy 
quarks through insertions of the SCET m Lagrangian. SCET operators that involve only light 
quarks, gluons, or at most one heavy quark correspond to standard fragmentation, whereas 
operators with a heavy quark-ant i- quark bilinear give rise to SD production. Since the fields 
in SCET have a power counting associated with them only a limited number of operators 
arise at each order in the power counting with corrections suppressed by A ~ wiq/pl- The 
fragmentation and SD production contributions are of different order in the SCET power 
counting, so we can consider each in turn. 

The generic form of the SCET operators that are required is 

O = O n O n O n *V^ Q O n , , (11) 

where 

^n',Q - \ Hn '>Q + Xn'){Hn',Q + X n '\ ■ (12) 

Here H n / Q is a quarkonium state that has a large light-cone momentum component Q in the 
n' direction, and X n > are states that are also collinear in the n' direction. The operators O n 
and On include fields in the n and n directions respectively that are collinear to the initial 
state, and O n > contains fields in the n' direction that will eventually hadronize into a jet 
that includes the quarkonium state. 

Given the four possible initial parton combinations (schematically qq, qG, Gq, GG), there 
are twelve types of fragmentation operators (fragmenting of a light quark, heavy quark, 
or gluon) and four types of SD operators. Since the steps we take to arrive at the final 

3 There is an additional scaling that we could consider. If Shad ^ -Epartonic then we are in the regime 
where the parton light-cone momentum fraction x goes to zero (a.k.a. small a;), and the existence of terms 
scaling like In x would be of concern. We will not consider such a scenario since it goes beyond the scope 
of this work. 
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FIG. 3: Matching Og„ at leading order. On the left is the square of the full theory amplitudes, and 
on the right is the matrix element of OrL At this order the matching coefficient is proportional to 
a 2 s {Q). 



factored form for the matrix element of each operator are the same, we will consider one 
operator of each type in detail. The results are easily generalized to the other situations. 
The fragmentation operator for an incoming qq to produce an outgoing gluon in the n' 
direction is 

0° = J dwicfc^C^, (13) 

where i, j, and k run from one to two. This operator is arrived at after Fierzing the full theory 
diagram. We have kept only the contributions that lead to a non-vanishing matrix element 
at leading power. This operator scales as A 6 and can produce a QQ pair through a time 
ordered product with an O(A ) interaction term from the SCET, m Lagrangian. The hard 
matching coefficient C(ui,Uj,u' k ) is determined by perturbatively matching this operator 
onto the full theory and is therefore given by an expansion in a s (Q). For example, the 
matching of O^L is depicted in Feynman diagrams in Fig. [3j The matching coefficient at tree 
level is proportional to a 2 (Q). However, since the production of a QQ requires the insertion 
of an interaction term from the SCET m Lagrangian, this operator will acquire an additional 
a s {2mQ). As a result the fragmentation contribution is proportional to a 2 s {Q)a s {2mo). 

The SD production operator that describes the production of a QQ pair from an initial 
light qq is 

°?? = J2 \ du; i dcu j duj' l C a (uj i ,u j ,u' l ) (14) 

X (Xn,o; 2 2Wi)(X<^r aM {l,T A }x<^^^ , 

where I = 1, . . . , 4, and T a ^ G \{f , ijl'j 5 , fY±>} with Y±_> = Y ~ n' v f/2 - n' u f/2. Note 
the heavy quark-ant i- quark pair can be in either a color-singlet or color-octet state. This 
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FIG. 4: Matching of Oqq onto the full theory at leading order. On the left are two leading order 
Feynman diagrams that contribute to the production of a QQ pair from an incoming qq pair . On 
the right is the tree- level matrix element of Oqq . The dashed lines are incoming and outgoing 
collinear light quark lines, and the dashed double lines are incoming and outgoing heavy quarks. 
At this order the matching coefficient is proportional to a 3 (Q). 



operator scales as A 8 in the SCET m power counting and is therefore A 2 ~ ^q/Vl suppressed 
relative to the fragmentation contribution. As for the fragmentation contribution, the hard 
matching coefficients C a are determined by perturbatively matching this operator onto the 
full theory. For example, the lowest order matching, shown in Fig. HJ is proportional to 

We pause here to discuss a point made in the introduction; namely that at leading order 
in a s the color-singlet 3 Si SD contribution scales as mg/p| relative to gluon fragmentation. 
This is in apparent contradiction with our operator analysis which leads us to conclude that 
the suppression of DPF is only of order A 2 ~ rnq/p 2 ^. The resolution is that at leading order 
in a s the color-singlet 3 Si SD contribution does not have a leading in A DPF component, in 
other words the matching of the color-singlet 3 Si SD contribution onto the leading DPFF 
vanishes at leading order in a s . However, there are power corrections to the DPFF that 
scale as A 4 (and higher) relative to the gluon fragmentation function. Such sub-leading 
DPFFs could for example have the form of the leading DPFF with factors of the SCET 
covariant derivative inserted or have explicit factors of the quarkonium mass. It is these 
subleading DPFF contributions onto which the color-singlet 3 Si SD contribution matches 
at leading order in a s . However, at NLO in a s , the color-singlet 3 Si SD contribution has a 
contribution whose scaling is vr?qjp\ suppressed relative to fragmentation and which matches 
onto a leading DPFF. 

Next we take matrix elements of the SCET operators between p states along the n and n 
directions. The usoft fields decouple from the collinear fields in the action by using the BPS 
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field redefinition [33] , which decouples the Hilbert spaces of the various modes and allows us 
to factorize the matrix element. After these steps the matrix element of the fragmentation 
operator in Eq. ( TT5|) is 

|0>. (15) 



x 



jjAfi ^BA^H V AC n C 
n n',oj' 1 I 1 n^Q 1 -°n» 



where Y AB are Wilson lines in adjoint representation that run along the light-cone from 
infinity to the point where the operator is situated, which we take to be the origin. Since 
the projection operator is limited to the collinear sector and the soft Wilson lines end at 
the same point as a consequence of the multipole expansion, the Y's cancel. Note that each 
state scales as A -1 and the projection operator V^, q scales as A~ 2 so the final matrix element 
of the fragmentation operator scales as A . As mentioned previously, it is easy to generalize 
this to the other initial and final states in SCET. 

Next we consider the matrix element of the SD production operator in Eq. ffbil) . After 
factoring usoft from collinear, this matrix element is 

(pnPn\0®Q\p n Pn) = / du i du) j du' l C a {uj i , U)j , u\) (p n | Xn,u 2 |xn, W1 \Pn) (Pn\ Xn^Xn^Pn) 

x (o|x„^^^^){l,^ A } Xn ,^p^Xn^^^^ ) {l,T A } Xn , M |o), (16) 

where again the soft Wilson lines cancel. The matrix element scales as A 2 so is suppressed 
relative to the fragmentation matrix element. 

The matrix elements involving the incoming states \p n ) and \p n ) in Eqs. ( TT5TT6]) are related 
to the parton distribution functions (PDFs) [34| : 

1 f 1 

^^(Pn{p))\Xn Ml iXn,u, 2 \Pn{p)) = 4fi • p dz 8{u -)8{u + - 2zfi ■ p) f i/p (z) (17) 

spin 

- 4n-p / dz5(u-)5(uj + + 2zn ■p)fi/ p {z), 
Jo 

l^tpnWMB^BW^ip)) = -^^J\z6(u-)8(u + -2zn-p)f g/p (z), 
spin 

where u± = u>i ±w 2 , fi/ p (z) is the quark PDF, fy p (z) is the anti-quark PDF, and f g / p {z) is 
the gluon PDF. Futhermore, the vacuum matrix element of the fragmentation operator in 
Eq. (fl5j) can be related to the standard fragmentation function that gives the probability of 
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finding in the gluon a quarkonium state H moving in the n' direction with large light-cone 
momentum n' ■ p: 

j^ T {o|'&K^A^J 1°) ( 18 ) 

c 

4 [ ldz xt >\z( 2fl '- p \ n ( \ 
= --T / —^-) S[z- —r- D H/g {z) . 

For completeness we give the SCET definition of the light-quark to quarkonium fragmenta- 
tion function, 

1 . , H _ . . f 1 dz , ( 2v! ■ p 



2N, 



M0\fXn>yV^ n ,. p Xn'y 2 \0) = 2 f -8{u'_)5(z - ^LP) D H/q (z) . (19) 

Jo z V UJ + J 



These definitions agree with those in Refs. 



23-1231 • Substituting Eqs. (fTTtlTS]) into Eq. (0SJ) 



we arrive at the familiar factored form for the fragmentation cross section in proton-proton 
collisions: 

d 2 a f dz 
^ 2 dp 2 dy = J dxidx2 ~^^ Xl ^ x ^ z ^^ y ^^/p^ x ^^/p^ DH /9^ ' ( 20 ) 

where a is the short-distance partonic differential cross section for producing a gluon from 
the collision of a quark and anti-quark. 

Double parton fragmentation is the kinematic situation in which a collinear, highly ener- 
getic, nearly on-shell heavy quark-anti-quark pair hadronizes into an energetic quarkonium. 
In every quarkonium production process, the heavy quark and heavy anti-quark have small 
relative momenta, so they are collinear to each other. What distinguishes double parton 
fragmentation from, for example, threshold production, is the large boost the heavy quark- 
anti-quark pair have relative to the lab frame. Because p± 3> mg we can think of the heavy 
quark and anti-quark as light-like collinear SCET modes. For this to have an invariant 
meaning, the heavy quark-anti-quark pair must be recoiling against one or more energetic 
jets of partons, so there exists a scale in the problem much larger than tuq. The energetic 
partons (quarks) will form a jet via collinear radiation and generate a set of large logs that 
would not be present if the quark pair was produced nearly at rest in the lab frame. 

Therefore, the vacuum matrix elements in the SD production operators are also fragmen- 
tation functions, but of a new type. The DPFF is defined in terms of the matrix element in 
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Eq. (USD by 



(OlXnX^Il.^^^pXnV^Il^l^lO) (21) 
= 8 -u'+uj'-uj')I — du dv S(z - f' P , )5(v - 1 - z^)5(u - z^-) 



X zD a{l,8}( U > V i Z )- 

This distribution is a combination of fragmentation function and light-cone distribution 
amplitude. The light-cone momentum fraction variables are 

n'-p n'-p ^ 



L0[ — U)' 2 Uj' 4 — L0' 3 

v = z—1- = 1 + z— 
n-p n-p 

u = z—^- = l + z—2-. 
n'-p n'-p 

The variable z corresponds to the fraction of the QQ pair light-cone momentum that H 
carries away. The variables u and v correspond to the fraction of the total QQ light-cone 
momentum carried by each of the heavy quarks in the QQ pair. These variables do not 
have to be the same. The only constraint on the momentum is that the difference of the 
total light-cone momentum of the two heavy quark-anti-quark pairs is zero. The expression 
above can be inverted: 

D®^ 8 y(u,v,z) = - j dco[duj2duj' 3 duj' 4 S(uj[ — co' 2 — — —)$( UJ 2 ~~ ~~ — {v — 1))5{uj' a — —^-u) 



z z z 



X 



(o|xn^^^ W{l,^ A }xn^^^n^ P Xn^a,i^^ ) {l,^ A }xn^^|o). (23) 



This definition of the DPFF is proportional to the one in Ref. [10|] with the following variable 
redefinition: u -> (1 + £)/2 and v -> (1 + C')/ 2 - Substituting Eqs. (H3 HQ into Eq. (ITo]) 
gives a generalized factored form for the SD production cross section 

d 2 a 1 f , , dz fn'-p\ 3 1 r lg . 
dydpl 2 J z \ z J (4vr) 2 

xf q/p (x 1 )fg /p (x 2 )D^ 18} (u, v, z) , 

where the factor of (n'-p/z) 3 comes from switching from the dimension-full variables Wj, u), 
and c^' to dimensionless variables X\, x 2 , z, u, v. Our result agrees with the factorization 
formula in Ref. \li 
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II. EVOLUTION EQUATIONS FOR THE DPFF 



In this section, we derive the evolution equations for the D^ 8 y(u, v , z) DPFFs. We 
consider these DPFFs in particular because they are most relevant to color-singlet 3 Si pro- 
duction. The diagrams for computing the one-loop anomalous dimensions are shown in 
Fig. [5J4.-G. In these diagrams the single lines represent SCET collinear fields and the double 
lines are Wilson lines. In addition to the seven diagrams, it is possible to generate additional 
diagrams by reflecting a diagram about the horizontal or vertical axes, or both. We refer 
to diagrams obtained from those in Fig. [5] by reflecting about the horizontal axis by adding 
a hat, e.g. B, diagrams obtained by reflecting about the vertical axis by adding a bar, e.g. 
B, and by doing both reflections by adding a hat and bar, e.g. B. Note that A = A, 
D = D, and E = E and these do not constitute distinct diagrams. The remaining diagrams 
have distinct images under the three possible reflections. As mentioned above we focus on 
those operators which have the Dirac structure ^'7^, since these are most relevant to the 
production of QQ pairs in 3 Si configurations. Note the this Dirac structure does not mix 
with other Dirac structures due to the symmetries of SCET. However, we allow for both 
octet and singlet operators since they do mix. 



A. Virtual Diagrams 

The divergent pieces of diagrams A and A vanish at leading power. If we consider the 
color-singlet operator and calculate diagrams B and C only, the IR divergences cancel and 
the UV divergence would lead to an anomalous dimension identical to the ERBEL 38 



evolution for light-cone wavef unctions. In general the infrared divergences in diagrams B 
and C do not cancel because the C diagram has a different color factor for the color-octet 
operator. Instead the IR divergences cancel between a number of different diagrams in a 
non-trivial way. 

In the diagrams in Fig. [5] the quark and anti-quark on the left-hand side of the cut have 
outgoing momentum and P3 respectively, while the quark and anti-quark on the right- 
hand side of the cut have incoming momentum p± and P2 respectively. We express the large 
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FIG. 5: The diagrams we need for the one-loop running. Not shown are the diagrams which 
are mirror images with respect to the horizontal and vertical axes. Diagram A reflected about a 
horizontal (vertical) is denoted in the text by A(A). 



components of these momenta in terms of momentum fractions: 



x 
A 



ri ■ (pi + p 2 ) 



(25) 



x- 



x- 



n ■ P4 

~P~ 

n' ■ pi 



where P is the large light-cone momentum component of the final state QQ pair. Another 
technical complication is that individual diagrams have rapidity divergences that cancel in 

the calculation. 

This is accomplished here by adopting the rapidity regulator of Refs. 
The result of evaluating diagram B is 

3 



es 


3 


30 


■ 



M B = 

Z7T e 



1 , zu 
- + ln — 



1 



J [jp ) 



where j = 1,8 indicates color-structure, and 



{l,T a Kn'^^{l,T a K, 



(27) 
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is a combination of SCET spinors in either a color-singlet or color-octet combination. The 
symmetric diagrams can be found via the simple replacements: 



M B = M B (u <->v,\<-> 0, 
M B = M B (u <-> u, A <-> A), 



= M B (u o v, A ^ f), 



(28) 



where w = 1 — w, t> = 1 — i>, A = 1 — A, and £ = 1 — £. We find that the sum is 



M s + M B + Mg + M 

a s C F 
2tv~T 



(29) 



\ uuvvP 11 



+ 4 



2P 



5(1 - */z)£(A - u)£(£ - w)O j . 



Now we consider diagram C and its reflections. These diagrams have different color 
factors depending upon whether or not the operator is color-singlet or color-octet. We 



denote the color factors by (3^' 8 \ where (3^ = Cp and (3^ 



2N C 



Diagram C yields 



Mi 



c 



2tt e 



1 ( zv 
- + ln — 

7] \wr 



(30) 



Again the diagrams related by symmetry can be obtained by making the substitutions in 
Eq. (12 8p . Note that diagrams B and C and their reflections do not lead to any mixing 
between singlet and octet operators. 

As before, the individual diagrams are not IR finite, and the result of summing diagram 
C and its reflections is 



M c + M c + M d + M d 
u9(u-X) 



4 

Z V 



- + In . 

7] \ UUVVP 4 " 



+ 



A(A-u)+ 



5{i - v) - 



2tt e 

v 0(f -v) v 9{v - 



5{\-u) 



6(X - u)6(Z - v) - 
2P J 



u6(\-u) 
A (A - u)+ 

5(1 - z/x)0 3 . (31) 



Notice that for the color-singlet operator, but not the color-octet operator, the rapidity 
divergences (i.e., the 1/r] poles) and the corresponding logarithmic terms cancel between 



diagrams B and C and their reflections. The remaining terms 
kernel that is similar to that of a light-cone wave function 38 



ead to an evolution equation 
391 ] . For the color-octet op- 



erator the rapidity divergences in diagrams B and C cancel against the rapidity divergences 
in the real emission graphs which we turn to now. 



1(3 



B. Real Radiation 



Now consider the real radiation coming from diagrams D — G and their reflections. We 
introduce the color factor matrices: 

Pn Pis 



P 



i.i 



a a I \ i Ni-2 ' ^ 1 i 

P81 P88 / \ 1 ~2N7 J \ 1 ~W C 



The first index, i — 1 or 8, refers to the color state of the initial and final state quarks in the 
diagram and the second index, j, refers to the color-structure of the operator. We present 
only the diagrams shown in Fig. [5j 



M 



D 



a s ( z \ 3 1 n x 2 l — z/x . / z -A r / z-. 

— — Pa^ — -^5 (v - -H 5 (u- -X) J 

2tt \2PJ e uv 3 z 2 A£ V x* J \ x 1 



a 



/ z \ 3 1 ■= x 2 1 — z/x 



« s / Z \ 3 1 

M F = — - ' I 



27r \2PJ e uv 



P 



5(1 - z/x) + 



1 /zi/ 
- + In — 
77 \uF 



ux 9(1 — z/x) 
Xz (1 — z/x)+ 



(32) 



M, 



G 



2?r V2Py et/y 



zv 



5(1 - z/x) + 



- + In 

X(J(v--^(J(tZ-^A) 



#(1 — z/x) 
Xz (1 — z/x) + 



The reflections of these diagrams can be determined by the same substitutions given in 
Eq. (1281) : As previously mentioned, the rapidity divergences cancel in the sum of the collinear 
diagrams. This is as it must be since the process we are considering is bereft of any soft 
sector which usually supplies the mechanism for the cancellation of rapidity divergences. 

III. RENORMALIZATION 



The DPFF is renormalized multiplicatively as follows 

f dz' 

D®(u,v,z) = / du'dv' ' —Zij(u,u' ,v,v' ,zj 'z', ii)Df(u\v' \z' pi). 

The renormalized distribution thus obeys 

d f dz' 

ji— Df (u" , v ", z", /i) = — / du'dv' — - r fij(u",u',v",v',z"/z',iJ,)Df(u',v',z',iJ,) 
dfj, J z 
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(33) 



(34) 



where the anomalous dimension is given by 



7i i (u,M / ,w,f / , z/z',/1) = (35) 
f dz" d 

J du"dv"—Zr 1 ( u , u", v, v", z/ z", fi)ii—Z aj (u" , u', v", v', z"/z', //) , 

and the indices % and j label the singlet = 1) and octet = 8) operators. The tree 
level matrix element of the DPFF using partonic states with momenta labelled by (x, A,£) 
is given by 



Dj(u, v, z) = (^) 3 5(1 - z/x)5(X - u)6(Z - v)O j 



(36) 



Given the results of the previous section and the wave function renormalization (which in 
SCET is identical to QCD), 

Z, = I - (37) 



we find the anomalous dimensions are given by 
7ll = -3^£Ls(u - u')5(v - v')5(l - z/z') 

7T 

a s C F „, z 



7T 



■6(v - v')5(l - -) 



z' 



6(v! — u)^-— * 



u 



v! [y! — u) + 



+ e(u-u')- 



a s C F 



S(u - u')5(l - -) 
n z' 



d{v'-v)- - 1 +Q(v-v'Y- 



u' (u — u') + 

V 1 



v' (v' — v) + 



v' (v — v')+ 



(38) 



781 = -V 9{1 - Z/Z,) ^) 2 \[(l-z/z% ' z/z' 



uv' + vu' 1 — z/z' 



—S(v-^v')5(u-^u') 



U'V 



+ 



uv' + vu' 1 — z/z' 



_(l-z/z') + z/z' 
uv' + vu' 1 — z/z' 

+ 



U'V 



_(l-z/z')+ z/z' 

uv' + vu' 1 — z/z' 
_(1 - z/z')+ + z/z' 



S(v - -v')S(u - - 



u 



—8{v - -v')5(u - -u') 

uv z z 

—S(v- - f v')8{u- -u') 
u'v' z' z' 



We also have 



718 



Cf 
2N, 



■7si 



(39) 
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Finally, 

788 = 



71 



a 1 



^-^-^(l--7) 



n2N L 
" vr 2N, 



■S(u - u')S(l - 4) 



9(u' — u) — -. : — h 9(u — u) — - 

J u'{u' -u)+ K J u'(u-u'). 

— V) — 7 r h V[V — V — 7 : — 

J v'(v'-v) + v J v'{v-v') + _ 



6{v 



a. 



-9(1 - z/z 1 ) 

71 



+ 



2N C 
uv' + vu' 
2N C [(I -z/z')+ 



uv' + vu' 1 — z/z' 

+ 



N?-2 



(1 - z/z')+ z/z' 
1 - z/z' 



U'V 



z/z' 



—5(v - -v')5(u - -u') 
uv z' z' 



1 

+ T7- 



uv' + vu' 



z/z' 



N c [(1 -z/z') + 
uv' + vu' 



1 



+ 



z/z' 
1 - z/z' 



-^8{v - - t mu - -,u>) 
u'v' z z 

±-,8{v - -v')5(u - - t u') \ . 

uv z z 



(40) 



—5{v--v')5{u--u) 



_(l-z/z')+ z/z' 

The anomalous dimensions computed in this section are new and are one of our main 
results. Evolving these equations from the scale p± to rriq will allow us to resum logs of 
P]_/mQ, but clearly solving these differential equations will be a complicated task that is 
beyond the scope of this paper. After evolving down to the scale rriQ, one must match the 
DPFF onto NRQCD matrix elements, which is the subject of the next section. 



IV. MATCHING ONTO NRQCD 

At the scale 2rriQ the heavy quark mass is integrated out by matching Df (u, v, z) defined 
in SCET m onto NRQCD. The SCET fields in the definition of this function contain light-like 
Wilson lines, which the NRQCD operators inherit. In the case where the NRQCD operators 
are in a color-singlet configuration the Wilson lines cancel and we arrive at a standard 
NRQCD long-distance matrix element. However, in the case where the NRQCD operators 
are in a color-octet configuration, the Wilson lines do not cancel. The presence of the Wilson 
lines in the color-octet operator ensures the proper infrared behavior of the operator matrix 



element 
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42]. However, because of the presence of the Wilson lines, these color-octet 
production operators are different from the ones introduced in Ref. 



To perform the matching of (u,v, z) onto NRQCD we adapt the framework for 



treating heavy quark effective theory in a boosted frame 



43|, 



4J to NRQCD in a boosted 



frame. Consider first a heavy quark field in NRQCD. In the heavy quark rest frame its 
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four-velocity is 

^ = (1,0,0,0) = ^ + ^. (41) 
The heavy quark momentum can be expressed as 

= m QV » + P + k", (42) 



where w M is the four- velocity of the heavy quark with v 2 = 1. Using the formalism of Ref. 45], 
the large components of the momentum become labels, denoted above by rnqv^ + k^, leaving 
derivatives acting on the field to scale as ~ rriQfi 2 , where the relative speed of the heavy 
quark and anti-quark (not to be confused with the heavy quark four- velocity v^) is /? <C 1. 
NRQCD gluons and massless quarks have momenta that scales as rriQf3 2 (i.e., they do not 
have labels). 

In a frame in which the heavy quark-anti-quark pair is boosted in the direction n' by a 
large factor Q, we have 

V* = IiL n 'M + I^'M (43) 

4m Q Q 

with a similar boosting for the other pieces of the momenta, such as 

Koost = \— n' ■ Kest ^ + ^n' ■ *v esi a'" + ^ , (44) 

where k rest are the components of momentum in the rest frame. So, for example, the NRQCD 
residual momentum scales as 



Q 

To match onto NRQCD we need to identify which components of the massive SCET mo- 
menta match onto components of the labels and residual momentum in NRQCD. A generic 
massive SCET momentum in the n' direction is p^ 1 = p M + r M , where the label momentum is 
ptj. = n' .pn'i 1 j2 + p^_, with scaling n'-p ~ Q, and pj_, ~ vtlq ~ QX. The residual momentum 
scales as r M ~ rriq/Q ~ QX 2 . Working in a frame where pj_, = 0, the label momentum 
n'-p in massive SCET matches onto the large component of the boosted velocity tuqu'-v 
in NRQCD. The SCET residual momentum will be split into label and residual NRQCD 
components. 

Next we consider how bilinears of massive SCET fields match onto bilinears in NRQCD. 
The generic SCET bilinear can be written as 

Xn',^ i{v) {XT A }Xn'^ = L>,p 2 W n ,6(n'-Vt - u 2 )r^{t,T A }6{n'-V - c^Wj^. (46) 
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When the virtuality drops below m 2 collinear parton splitting is no longer possible. However, 
the collinear Wilson lines in SCET m match onto another set of Wilson lines which can 
be thought of as arising from boosting ultra-soft Wilson lines that arise in NRQCD after 
performing a BPS field redefinition. We call such gluons ultra-collinear and designate the 
field as A£ c (x). Note that the momentum scaling of these gluons is down by a factor of ft 2 
compared to the usual collinear SCET gluons. This is because in matching we have reduced 
the virtuality of the external states from m 2 to m 2 /3 4 . As a result we find that the SCET 
Wilson line matches onto an NRQCD Wilson line 



1 9 



ds n'-A uc (n's) 



W n >(x) W uc {x) = Pexp 
In the matching the delta functions in Eq. (146 p fix u: 

8{n'-V-u)Wl,{x)^ 3 {x) 8{n'-V-u)Wl{x)^ 3 {x) 

= 5{m Q n'-v -u)Wl c {x)£ n > 3 {x) . 



(47) 



(48) 



We were able to push the delta function past the ultra-collinear Wilson line since n'-V acting 
on the Wilson line is again down by f3 2 compared to when it acts on the quark field. 

Finally we consider the matching of the bilinears of quark and anti-quark fields. For this 
step we can ignore the Wilson lines and delta functions, and focus on 



^ l2 rM{i,T A }^, Pl . 



(49) 



We reintroduce the large phase and relate the SCET fields to standard 4-component QCD 



fields 



3l|: 



-ip-x 



£n',p(-£) 



in' -d 2 J 



2p(x). 



Inserting this identity into Eq. ( 1491) and simplifying we find 



-i(p 1 -p 2 )-x 



^^{1,^}^, = ^x)r^{i,T A }^x) . 



(50) 



(51) 



Matrix elements of the bilinear in terms of QCD fields can then be related to matrix el- 
ements of bilinears of NRQCD fields using the techniques of Ref. 46|], with the matching 
schematically given by 



tP(x)r^{l,T A }^(x) ^ Ct\i^(a,D){l,T A }^ v + h.c. 



(52) 
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where C| are matching coefficients, ip v and \ v are the heavy quark and anti-quark fields 
respectively, a are the Pauli matrices, D is the covariant derivative in NRQCD, and £(cr, D) 
is a fixed function of a and Data given order in the NRQCD expansion. We have suppressed 
vector indices on and D). Putting all these pieces together, the matching from 
SCET onto NRQCD has the form 

X n .^ l{u) {^ TA }Xn>^ (53) 
-)• c\ v) 8{m Q n'-v - u x )5{m Q n'-v + w 2 )xiW MC E(a,D){l,T A }W^ c ^ + h.c. 

For a color-singlet configuration W uc Wl c = 1 and the Wilson lines cancel in NRQCD. 
We now match a generic DPFF onto NRQCD operators 

£(8) («,«,*) (54) 

„ r 2mnn'-v\ ( _. 2mnn'-v, ,\ / _. 2mnn'-v 
— >• Pj2; o I Zirtqu -v I o I mgn H (t> — 1)1 o I -U w 



z J \ z J \ z 

x (0| X ^ uc E(a, D){1, ^Wt^-VlWucS^, T A }^ cX „|0) 

x (0|x^ c ^D){l,T A }J^ , 
where X>j are matching coefficients, and 

V H " = J2\Hv + X uc ) (H v + X uc \ , (55) 

with if„ the NRQCD quarkonium state and X uc ultra-collinear states. As pointed out above, 
if the heavy quark bilinear is in a color-singlet configuration all ultra-collinear Wilson lines 
cancel and we are left with the standard NRQCD production matrix elements. However, if 
the heavy quark bilinear is in a color-octet configuration, the ultra-collinear Wilson lines do 



not cancel and we obtain co 
as first proposed in Refs. 41 



or-octet NRQCD matrix elements with light-like Wilson lines 



42]. 



As a concrete example let us consider the matching of D®®(u,v,z) to the color-singlet 
3 Si operator at leading order in the NRQCD expansion. The Dirac structure we need to 
match is 



Q 2 

— 5(m Q n'-v - u)[)5{mQn' -v + ui' 2 )5(m Q n' -v + Lo' 3 )6(m Q n' -v - u' i )xl(T l ^ v T >Hv il)la l x v . 

LlrriQ 
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Thus, 



D? Q (u,v,z) -J—S(l-z)5(l-v)5(l-u)(0\xla^ v V H ^yx v \0). (57) 



YlQrriQ y ' \2 J \2 

If running is neglected and Eq. (1ST)) is inserted into Eq. ( 124)) then we recover the results of 
the NRQCD factorization formalism. 



V. CONCLUSIONS 



In this paper we used SCET to derive a factorization formula for fragmentation and double 
fragmentation production of quarkonium. In addition, we have shown that production in 
the regime where the hard scattering scale is much larger than pj_ is suppressed. We have 
presented for the first time the anomalous dimensions of the two by two system of DPFFs 
that are relevant for the production of 3 Si states and showed how to match the DPFF onto 
NRQCD production operators. The running of the DPFF sums logs of p±/mq and could 
have a large effect on the quarkonium production rate. A future publication will calculate 
the effects of the running of the DPFF on both the production rate and the polarization. 
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Appendix A: SCET review 



SCET is an effective field theory coupling soft and collinear degrees of freedom. Collinear 
degrees of freedom have light-cone momenta (k + , k~, k±) that scale as 

p~Q(A 2 ,l,A), (Al) 

where Q is the large mass scale and A is the SCET expansion parameter. Soft modes have 
momenta that scale as 

P~Q(A,A,A), (A2) 

while ultra-soft (usoft) modes scale as 

p~Q(A 2 ,A 2 ,A 2 ). (A3) 

For example, if we are interested in describing the motion of a highly energetic particle 
with off-shellness m 2 in the light-like direction n, the collinear mode light-cone momentum 
will scale as Q(X 2 , 1, A), where A ~ m/Q with m <^ Q and the usoft mode momenta scale 
as (A 2 , A 2 , A 2 ). Which modes are present depends on the process. In order to keep the 
parametrically different momenta separate we introduce a projection operator V 11 which 
projects out momentum of order Q or QX. The derivative only operates on residual 
momenta ~ QX 2 . 

SCET operators are constructed out of gauge invariant combinations of fields and collinear 



light-like Wilson lines 34j which for convenience are combined into a single field. For quarks 



that have large momentum component in the n direction we define 

Xn ^= [S(u-V)W^], (A4) 

where the quark field £ n) p is labelled by the light-cone direction n and by the large light-cone 
momentum components of p: p = n-pn^/2 + pj_. The operator V = n ■ V acts on the 
quark field to project out the large light-cone momentum label: V^ n ^ = n-p^ n ,p- The SCET 
collinear Wilson line is defined as 



Wl 



(A5) 



. perms 

where A% is the collinear gluon field. The gauge invariant field strength is 



! 



W = -- [6(u - V)Wi\iV» + gA^ q , iVl + gAl Q ]W n \ (A6) 
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where 

— fi^ 1 

iK = y v + V ± + Y n ' {w + gAus) ' (A7) 

with A us being the ultra-soft gluon field. Since (Gn,u) tl ' U is not homogeneous in the SCET 
power counting we project out the leading contribution by introducing the gauge invari- 
ant field VB% U = in u (Q nLU ^ ,L '<7_Lu; "where g^ = g^ u — n^n u /2 — n fl n u /2 projects out the 
components of a four vector that are perpendicular to n and n. Collinear operators are con- 
structed out of combinations of the gauge invariant fields above. For example, an operator 
that creates a quark and an anti-quark moving in opposite light-like directions is Xn,Q^Xn,ui, 
where T is a direct product of a Dirac matrix, and a color matrix depending on the current 
producing the quark-anti-quark pair. The form of V is constrained by the symmetries of 
SCET and for this combination of fields is restricted to T = {I, T a } <g> {1, 75, 7^}, where 
7^ = 7i,<7j^ '. Since the production of two back-to-back light-like quarks with large energy is 
associated with a current that has large invariant mass the operator in our example must 
be matched onto QCD. This gives the matching coefficient which can depend on the large 
light-cone momentum components u and u. As a result the short- distance coefficient and 
operator are convoluted : 

QCD -> J (kudu C(u, Ld)Xn,u>FXn,uj ■ (A8) 

Furthermore, usoft modes can be decoupled from collinear modes in the SCET Lagrangian 
through a field redefinition 33 ] 

n,q T ± n- r> -n,q ± n 

W n -»• Y n W n Yl 

Xn,u) ^ Y n Xn,u) 1 (A-9) 

where the collinear fields on the right no longer couple to usoft fields, and Y n is a path-ordered 
exponential of the usoft gluon field 

Y n (x) = P exp^ig J ds n- A us (sn + x)^j . (AlO) 
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